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1. Introduction. Reduction problem for Feynman integrals 

After a tensor reduction is performed a given Feynman graph F results into various scalar 
Feynman integrals that have the same structure of the integrand with various distributions of powers 
of propagators 

f [d^kx..A'^kh 
F{ai,...,an) = J ■■■ J (1.1) 

d = 4 — 2£; the denominators are either quadratic or linear with respect to the loop momenta 
Pi = ki, i = I,... ,h and to the independent external momenta p^+i = qi,... ,ph+N = of the 
graph. 

There are different methods of evaluating those integrals. An old analytical strategy is to 
evaluate, by some methods, every scalar Feynman integral generated by the given graph. And 
what has already become a traditional strategy is to derive, without calculation, and then apply 
integration by parts (IBP) relations between the given family ofFeynman integrals as recurrence 
relations. 



After such a reduction a general integral from the family ( |1.1D is expressed as a linear combi- 
nation of some basic (master) integrals. Thus, the whole problem of evaluation falls apart into two 
parts: constructing a reduction procedure and evaluating master integrals. This paper is devoted to 
the first part of the problem. 

Let us introduce some notation. Let ^ = F{ai, . . . ,an) be functions of integer variables 
a 1, it is an infinitely dimensional linear space. The simplest basis consists of elements 

Hai,...,a„ia'l,---,a'n) = Sai,a[ ■■■^a„,a'„- (1-2) 

From this point of view, Feynman integrals form a point in all relations we use turn into 
linear functional on this vector space. 

There are different types of relations. The most commonly used among them are IBP relations 

101 

/.../d%d^fc...jL(p,j^)=0, (1.3) 

After differentiation we obtain some relations of the following form: 

Y^aiF{ai+bi^u...,an + bi^„)=0. (1.4) 

There are more relations one can consider: Lorentz-invariance (LI) identities [^], symmetry 
relations, e.g., 

. ..,«„) = (-l)'^>'"+-''"''"F(a^(i), . . . (1.5) 
boundary conditions: i.e. the conditions of the following form: 

F{ai,a2,. ■■,an) =0 when a,, < 0, . . .a,j < (1.6) 
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for some subsets of indices ij; parity conditions and others. All those relations lead to a possibility 
to express given Feynman integrals in terms of other Feynman integrals. 

Therefore one has to name certain integrals irreducible (master) and aim to reduce any other 
integral to those. An attempt to formalize the definition of master integrals was made in [p^. To 
define master integrals one needs to introduce a priority between the points (ai, ... ,£?„), hence an 
ordering. 

There are different ways to choose an ordering. To start with, Feynman integrals are simpler, 
from the analytic point of view, if they have more non-positive indices. Moreover, solving IBP 
relations by hand usually consists of reducing indices to zero. 

Naturally, we come to the notion of sectors. Those are 2" regions labeled by subsets v C 
{1, . . . ,n}, where Oy = {(ai, . . . ,a„) : c?,- > if / E V , a, < if / ^ v}. A sector Ov is said to 
be lower than a sector if V C jJ.. Furthermore, , . . . >- F(aj , . . . if the sector of 
{a\,...,a'^) is lower than the sector of (ai , . . . , To define an ordering completely one has to 
introduce it in some way inside the sectors (this will be discussed below). 

Initially they used to solve relations by hand, but with the growth of the complexity of the 
problems it turned to be impractical. The first algorithmic approach to solving IBP relations was 
the so-called Laporta's algorithm [^,[TC, 13, 14] which is now very well known and actively used in 
practice. It is based on the fact that the total number of IBP equations written for concrete indices 
grows faster than the number of independent Feynman integrals, hence this system of equations 
sooner or later becomes overdetermined, and one obtains a possibility to perform a reduction to 
master integrals. At the moment, there is one public implementation of this algorithm called AIR 
\^ and a lot of private versions. ' 

There are other algorithmic approaches to solving the relations such as the Baikov's method 
^, ||] and the use of Grobner bases The reduction using Grobner bases has been suggested 



by Tarasov pO| , 2 1 1 who reduced the reduction problem to differential equations by introducing 
a mass for every line, a,i^ A direct application of Grobner bases (without the use of 

differential equations) was initially suggested by Gerdt [11]. 

One more approach initially based on Grobner bases, called the 5^-bases, has been developed 
in [^] (see also [ p^ for a brief review). It was applied in practice for a reduction of a family of 
three-loop Feynman integrals necessary for the analysis of decoupling of c-quark loops in Z>-quark 
HQET [12]. In the next section we outline basic ideas of this approach in an updated form. Then, 
in Section 3, we formulate how it can be combined with other strategies, in particular with the 
Laporta's algorithm. 



2. The 5-bases approach 

Suppose first that we are interested in expressing any integral in the positive sector aji „| as 



a linear combination of a finite number of integrals in it. The left-hand sides of IBP relations (1.4) 
can be expressed in terms of operators of multiplication A,- and shift operators = i+ , Yf = i , 



^by T. Gehrmann and E. Remiddi, M. Czakon, P. Marquard and D. Seidel, Y. Schroder, C. Sturm, A. Onishchenko, 
O. Veretin, . . . 
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where 

(A-F)(ai,...,a„) = a,-F(ai , . . . and (2.1) 
{Yi^ ■F){ai,...,a„) = F(fl:i, . . . ,a,- ± 1, . . . ,£?„) . (2.2) 

Let M,...,„ be the algebra generated by shift operators Y^^ and multiplication operators A,. 
It acts on the field of functions ^ of « integer variables. One can choose certain elements fi 
corresponding to IBP relations and write 

{frF){ai,...,a„)=0. (2.3) 

The choice of elements fi is not unique, we will choose them so that they do not depend on 
Yi^. Consider the left ideal C £/i n, generated by the elements /,-. This ideal is named as the 
ideal of IBP relations. For any element X G we have 

(XF)(1,1,...,1)=0. (2.4) 

Also we have 

Fiaua2,...,an) = {Y,"'-' ...Y^"-'F)il,l,...,\). (2.5) 



The idea of the algorithm is to reduce the operator on the right-hand side of (2.5) using the 
elements of the ideal J^. Suppose we are interested in F(ai,a2, . . . The reduction problem 
becomes equivalent to reducing the monomial F"'^' . ..Y""^^ modulo the ideal of the IBP relations. 
After obtaining an expression like 

Ff . ..Y,r' = In/; + Ic,,...,,Fr > . ..Yt' (2.6) 



it is left to apply ( |2.6| ) to F at ai = 1 ,...,«:„ = 1 and obtain the following expression: 

F(fli,a2, . . . ,a„) = Ic,-,,...^;,F(/i,/2, . . . ,/„) , (2.7) 



where integrals on the right-hand side are "master integrals" (the formulas (2.3) and (2.5) are used). 

To do the reduction we need an ordering of monomials of operators F, or, similarly, an or- 
dering of points (ai,...,a„) in the sector: (for two monomials Mi = Fj'^' . . .F„'"^' and M2 = 
F/'"^..F/""' we have (Mi • F)(l, . . . , 1) y (M2 •F)(l, . . . , 1) if and only if Mi ^ M2). Then the 
reduction procedure becomes similar to the division of polynomials. But one needs to introduce a 
proper ordering. 

A monomial is defined by its degree, i.e. a set of n non-negative integers (N"). Thus defining 
an ordering on monomials is equivalent to defining an ordering on N". We require the following 
properties: 

i) for any a € N" not equal to (0, . . . 0) one has (0, ... 0) -< a; 

ii) for any a,b,c gN" one has a ~<bif and only if a + c ~< b + c. 

For example, a lexicographical ordering is defined the following way: a set (/i, . . . ,/„) is said 
to be higher than a set {j\ (denoted by (/i , . . . , /„) (71 , • • • ,in)) if there is / < « such that 
'1 = 71, h = ii,..., ii-i = ji-i and // > j,. 
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Another example is a degree-lexicographical ordering: (/i , . . . , /„ ) y {ji,...,j„) if >LA, 
orL4 = LA and (/i,. ..,/„) >- (71,..., 7«) in the sense of the lexicographical ordering. 

An ordering can be defined by a non-degenerate n xn matrix {ak.i)'- for two sets of numbers 
(/i , . . . , /„) and (ji , . . . , jn) one first compares Y,i k<^\.i and ii<^\.i- If the first number is greater, 
then the first degree is greater; if the first number is smaller, then the first degree is smaller; and if 
those numbers are equal we compare Y,i Uai,! and Y,! ji'^ij and so on. 

The following matrices correspond to a lexicographical, a degree-lexicographical and a reverse 
degree-lexicographical ordering for n = 5: 

/l 0\ /l 1 1 1 l\ /l 1 1 1 l\ 

01000 10000 11110 

00100 , 01000 , 11100 

00010 00100 11000 

Vooool/ yoooio/ \iooooy 

Such an approach encounters the following problem: the reduction does not always lead to 
a reasonable number of irreducible integrals, so one has to build a special basis of the ideal first. 
Obviously having elements with lowest possible degrees is equivalent to obtaining master integrals 
with minimal possible degrees. Therefore one needs to build special bases. This can be done by an 
algorithm based on the Buchberger algorithm with the use of S'-polynomials and reductions [^. 

Moreover, one must keep in mind that we are interested in integrals not only in the positive 
sector. Our algorithm [ [T^ aims to build a set of special bases of the ideal of IBP relations (5-bases). 
The basic idea is to consider consider the algebra £/v generated by operators A,- and operators Y^^ 
for / E V and Yf' for / v. Then for Oy we again consider the ideal of IBP relations in £/v Now one 
has to construct sector bases (5-bases), rather than Grobner bases for all sectors, where an 5-basis 
for a sector Oy is a set of elements of which provides the possibility of a reduction to master 
integrals in Oy and integrals whose indices lie in lower sectors, i.e. Oy' for v' C V. 

This leads to considering many sectors — seemingly the problem becomes harder. But the 
important simplification is that one is not trying to solve the reduction problem in each sector 
separately but allows to reduce the integrals in a given sector to lower sectors — similarly to the 
"by hand" solutions. It is also worth noting that it is most complicated to construct s-bases for 
minimal sectors. 

The construction of 5-bases is close to the Buchberger algorithm but it can be terminated when 
the "current" basis already provides us the needed reduction. The basic operations are the same, 
i.e. calculating 5-polynomials and reducing them modulo current basis, with a chosen ordering. 

After constructing 5-bases for all non-trivial sectors one obtains a recursive (with respect to the 
sectors) procedure to evaluate , . . . ,a„) at any point and thereby reduce a given integral to mas- 
ter integrals. A description of the old version of the algorithm (implemented in Mathematica) 



can be found in [ 15]. 



3. Combining 5-bases with other strategies 



Constructing a basis requires a certain skill to find a proper ordering. At a high level (starting 
from 10 to 12 indices) it turns to be a certain problem. In some sectors the ordering can be found 
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automatically, in some sectors it is easily done manually, but in some sectors the 5-bases can be 
hardly constructed. There are several ideas that might allow the algorithm to construct the bases 
better and we hope that they will be implemented in the future. Currently this is not a problem of 
cpu time or memory, but an algorithmical one: with a fixed ordering the algorithm either produces 
a basis quickly, or does not do it at all. 

Recently we have worked out one more way to improve the algorithm. Experience shows 
that j-bases are constructed easily if the number of non-positive indices in a given sector is small. 
And what can help is the fact that if the number of non-positive indices is large, there is usually a 
possibility to perform integration over some loop momentum explicitly in terms of gamma func- 
tions for general d. However, to derive the corresponding explicit formulae, with multiple finite 
summations, for all necessary cases turns out to be impractical. 

In this situation, there is another alternative. Let us distinguish the propagators (and irreducible 
numerators) involved in such an explicit integration formula. They are associated with a subgraph 
of the given graph. Let us solve IBP relations for the corresponding subintegral in order to express 
any such subintegral in term of some master integrals. 

Then, after using this reduction procedure, it will be sufficient to use explicit integration for- 
mulae only for some boundary values of the indices. This replacement is very simple, without 
multiple summations. Coefficients in these reduction formulae depend not only on d and given 
external parameters but also on the propagators of the given graph, which are external for the sub- 
graph. 

It turned out possible to implement the solution of the recursive problem for the subgraph in 
terms of the Feynman integrals for the whole graph. In this reduction, pure powers of the parame- 
ters which are external for the subgraph transform naturally into the corresponding shift operators 
and their inverse. Integrals which are obtained from initial integrals by an explicit integration over 
a loop momentum in terms of gamma function usually involve a propagator with an analytic regu- 
larization by an amount proportional to e. Our code is applicable also in such situations. 

Even if we have not constructed the j-bases for all sectors we can still run the reduction pro- 
cedure. However, in some sectors there will be an infinite number of irreducible integrals left. 

One of the ideas is to use the Laporta's algorithm in those sectors. Still it leads to several 
problems in linking two algorithms together. One of the disadvantages is that the reduction works 
much slower if the number of irreducible integrals goes high. 

Another idea is to combine two algorithms inside one code. The 5-bases might give the reduc- 
tions in all sectors, where a basis can be constructed. And for the remaining sectors the systems of 
initial EBP's are used like in the Laporta's algorithm. 

Let us summarize the ideas of this paper. The j-bases approach seems to be a good addition to 
the methods of automatic integral reduction. The bases construction part of the algorithm should 
still be improved, but even now it is capable of constructing bases for quite complicated cases. A 
constructed basis in a sector means that in this sector one needs no more real solving of equations 
— the system of linear equations becomes upper-triangle. 

Moreover, it looks like the future is in combining the j-bases and the Laporta's algorithm. An 
algorithm named FIRE (Feynman Integral REduction) combining all mentioned ideas is currently 
being developed. We plan on performing much more complicated three-loop calculations with the 



6 



Reduction ofFeynman integrals 



A.V. Smimov 



use of FIRE (together with M. Steinhauser; the project is in progress). Afterwards we intend to 
make the algorithm pubUc. 
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through grant 05-02-17645. 
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